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In 1985, Beeman, Tsu and Thorpe established an almost 
hnear relation between the Raman transverse-optic (TO) 
peak width T and the spread in mean bond angle A8 in a-Si. 
This relation is often used to estimate the latter quantity in 
experiments. In the last decade, there has been significant 
progress in the computer generation of sample networks of 
amorphous silicon. Exploiting this progress, this manuscript 
presents a more accurate determination of the relation be- 
tween T and A8 using 1000-atom configurations. Also inves- 
tigated and quantified are the relations between the TO peak 
frequency and the ratio of the intensities of the transverse- 
acoustic (TA) and TO peak, both as functions of A9. As 
AO decreases, the TA/TO intensity ratio decreases and the 
TO peak frequency increases. These relations offer additional 
ways to obtain structural information on a-Si from Raman 
measurements. 

I. INTRODUCTION 

Many structural properties of a-Si, such as defect con- 
centration and variation in mean bond angle, are difficult 
to determine experimentally. This is because it is impos- 
sible to measure directly the coordinates of the atoms 
in a-Si. However, important information on the struc- 
ture of a-Si can be obtained indirectly, through a num- 
ber of experimental techniques. These techniques include 
neutron, x-ray and Raman scattering, electron-spin res- 
onance, and x-ray photo-absorption. Compared to other 
methods, Raman scattering is more sensitive to small 
changes in the short-range order of a-Si. For this rea- 
son, Raman measurements on a-Si are frequently used to 
obtain structural information [1-4]. 

The experimental Raman spectra of a-Si show two dis- 
tinct peaks, at about 150 cm -1 and 480 cm -1 , associated 
with the transverse acoustic (TA) and the transverse op- 
tic (TO) vibrational modes, respectively. Certain fea- 
tures in the Raman spectrum are highly sensitive to the 
structural properties of the a-Si sample. For example, the 
width of the TO peak is related to the root-mean-square 
bond-angle variation A9 in the amorphous network [5] . 

In several computational studies [6-8] , the relation be- 
tween r and AO was quantified. All studies indicate a 
broadening of the TO peak with increasing AO, but there 



is no quantitative agreement. Bccman's linear relation, 
r = 15 + 6A0, which dates back to 1985, is often used by 
experimentalists to determine AO from Raman measure- 
ments. Here, V is in cm -1 and AO in degrees. 

Beeman derived his relation using nine structural mod- 
els of a-Si. Of these models, five were generated from 
the same 238-atom, hand-built model by Connell and 
Temkin [9], which contains only even-membered rings. 
In contrast, all simulations on a-Si find an abundance of 
five- and seven- fold rings. Moreover, these five Connell- 
Temkin models are statistically dependent and not pe- 
riodic, consequently containing a large fraction of sur- 
face atoms. Experimental values of AO, based on the 
radial distribution function of a-Si obtained in neutron- 
diffraction studies, range from 9.9 to 11.0 degrees [10]. Of 
the nine structural models used by Beeman, only three 
exhibit values of AO in this range. New techniques to 
generate a-Si structures, such as ART [11-13], as well as 
more powerful computers, have made it possible to gen- 
erate larger and more realistic a-Si systems via computer 
simulation. 

Also the description of the Raman scattering pro- 
cess has improved. Beeman used the bond polarizibility 
model proposed by Alben et al. [14], which dates back 
to 1975. Characteristic for this model is the inclusion 
of three weighting parameters, whose values must be set 
somewhat arbitrarily. Several studies have indicated that 
the values originally proposed by Alben yield an incor- 
rect value for the depolarization ratio [3,4]. These stud- 
ies therefore propose different weights. Since then, other 
polarizibility models have been proposed, for example by 
Marinov and Zotov [15]. 

In this manuscript, we will re-investigate the relation 
between T and AO by computer simulation. This simu- 
lation is based on a large number of 1000-atom, periodic 
configurations, with structural properties (radial distri- 
bution function, spread in mean bond angle) that are 
in excellent agreement with experiment. Furthermore, 
recent advances in neutron scattering techniques have 
made it possible to directly compare the bond polariz- 
ibility models to experiment [16]. We therefore also in- 
clude a detailed comparison of the model of Alben and 
the model of Marinov and Zotov to experiment. 
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Additionally, we present two other methods to obtain 
structural information from the Raman spectrum. The 
TA/TO intensity ratio [1] and the location of the TO- 
peak [17] are believed to be directly related to A9; these 
relations will be quantified. 

The outline of this paper is as follows. In section II, 
we explain the generation of the a-Si configurations used 
in this study. We then discuss how the Raman spectrum 
is obtained from these configurations. The results and 
conclusions are presented in sections III and IV, respec- 
tively. 



3. Every 50 ART moves, up to approximately five 
ART moves per atom when the energy has reached 
a plateau, the configuration is stored. For 1000- 
atom samples, this procedure yields a set contain- 
ing 100 samples. 

This procedure is repeated nine times, generating nine 
statistically independent sets of 100 correlated configura- 
tions each. For each set, we found that A6* ranges from 
approximately 10° for the well-annealed configurations, 
to approximately 14° for the poorly annealed configura- 
tions. 



II. METHOD 

To calculate Raman spectra, three ingredients are re- 
quired: (1) a potential describing the atomic interactions 
in the sample, (2) a continuous random network repre- 
senting a realistic sample of a-Si, and (3) a model assign- 
ing Raman activities to the vibrational cigenmodes of the 
sample. 

In the present work, we use a modified version of the 
Stillinger- Weber potential for all calculations. This po- 
tential has the same functional form as the original SW 
potential [18], but with different parameters. The pa- 
rameters were chosen specifically to describe a-Si, see 
Rcf. [19]. 



A. Sample generation 

To study the effect of A8 on the Raman spectrum, we 
require a number of a-Si configurations with varying val- 
ues of A9. These configurations were generated using the 
activation-relaxation technique (ART) [11-13]. As was 
shown in previous studies [11,12], ART yields structures 
in good agreement with experiment. They display a low 
density of coordination defects, a narrow bond-angle dis- 
tribution and an excellent overlap with the experimental 
radial distribution function (RDF). The method is out- 
lined below: 

1. Initially, 1000 atoms are placed at random in a pe- 
riodic cubic cell; the configuration is then relaxed 
at zero pressure. 

2. The configuration is annealed using ART. One 
ART move consists of two steps: (1) the sam- 
ple is brought from a local energy minimum to a 
nearby saddle-point (activation), and (2) then re- 
laxed to a new minimum with a local energy mini- 
mization scheme including volume optimization, at 
zero pressure. The new minimum is accepted with 
a Metropolis probability at temperature T = 0.25 
eV. 



B. Calculation of Raman spectra 

We focus on the reduced Raman spectrum I(u), with 
thermal and harmonic oscillator factors removed. This 
spectrum is a function of frequency u> and of the form 



/(w) = C(uj)g(uj), 



(1) 



where g(u>) is the vibrational density of states (VDOS) 
and C{ui) a coupling parameter, which depends on fre- 
quency and on the polarization (HH or HV) of the inci- 
dent light used in the Raman experiment. 

To calculate the VDOS, the hessian is calculated. Di- 
agonalization of the hessian gives the frequencies of the 
vibrational modes, from which the VDOS is obtained. 

The function C(u>) is obtained from the polarizibility 
tensor a(ui). In terms of a(uj), the coupling parameter 
for HH and HV Raman scattering becomes Chh(u p ) = 
7G 2 + 45 A 2 and C H v{u p ) = 6G 2 , respectively, with A 
and G 2 given by the tensor invariants 



A = - [an + a 22 + a 33 ] 



and 
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see for instance Ref. [20]. 

The form of the polarizibility tensor a(ui) still needs 
to be specified; this is the most uncertain part of the 
calculation. Several models have been proposed, amongst 
which the commonly used model of Alben et al. [14] and 
the more recent model of Marinov and Zotov [15]. 

In the model of Alben, a cylindrical symmetry of the 
individual bonds is assumed and each bond is treated in- 
dependently as a homopolar, diatomic molecule. Three 
different forms for the bond polarizibility tensor are in- 
troduced: 
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Here, the summation runs over all atoms I in the sample 
and their nearest neighbors A, r\ is the unit vector from 
the equilibrium position of atom / to the nearest neigh- 
bor A, ui is the displacement vector of atom I when it 
is vibrating in mode p and I is the unit dyadic. The to- 
tal polarizibility tensor a is a weighted sum of the three 
terms, i.e. a = B\ol\ + B2CX2 + B^a^. As was stated in 
the introduction, the precise choice of the weights Bi is 
somewhat arbitrary and this is the major shortcoming of 
the model. Several studies have indicated that mecha- 
nisms 1 and 3 provide the main contribution to the Ra- 
man scattering process; these propose to use B\ : B2 : B3 
proportional to 2 : : 1, respectively [3,4]. In this paper, 
we will use this set of weights. 

The model of Marinov and Zotov (MZ) has no free 
parameters. In this model, the bond polarizibility is 
expressed as a sum of three components; two compo- 
nents parallel to the bond arising from bonding and non- 
bonding electrons and a third component perpendicular 
to the bond, see Ref [15]. Under these assumptions, the 
polarizibility tensor takes the form: 



+ ^ (b m f m +r m b. 
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Here, the summation runs over all bonds m in the sample, 
f m is a unit vector parallel to the bond, r m is the bond 
length and b m is defined as Uj — uf, where Ui and Uj are 
the displacement vectors of atoms i and j constituting 
the mth bond, when vibrating in mode p. 

The coupling parameter for a-Si has also been de- 
termined experimentally [16] through neutron scattering 
methods. According to this experiment, the coupling pa- 
rameter is a slowly increasing function of frequency. We 
will use the experimental result of Ref. [16] to test the 
validity of both the model of Albcn and the MZ model. 



III. RESULTS 

First, in section III A, we compare the polarizibility 
models of Alben and Marinov and Zotov to experiment. 
In the subsequent sections, we investigate the relation 
between the spread in the bond angle AO and: 

1. Raman TO peak width, 



2. Raman TO peak position, and 

3. Raman TO/TA intensity ratio. 



We show results for HV polarized light only; this is the 
usual experimental situation. Results for HH polarized 
light have also been obtained and are available upon re- 
quest. 



A. Raman coupling parameter 

The solid curves in Fig. 1 show the HV Raman coupling 
parameter for a-Si calculated using the model of Alben 
(top) and the MZ model (bottom) for the bond polariz- 
ibility. The experimental result of Ref. [16] is also shown 
(dashed). For this calculation, we used a well-annealed, 
1000-atom configuration with AO = 10.0°, since this will 
most closely resemble the experimental sample. We have 
checked that the general features of the curves in Fig. 1 
do not depend on the details of the configuration used: a 
number of other, well-annealed, configurations, with A9 
ranging from 10.0° to 11.0°, gave similar results. 




100 200 300 400 500 600 




100 200 300 400 500 600 



FIG. 1. HV Raman Coupling parameter for a-Si calculated 
using the model of Alben model (top) and the MZ model 
(bottom). The dashed line is the experimental result taken 
from Ref. [16]. Frequency is in cm -1 and all curves are 
area-normalized. 

Fig. 1 shows that both models yield an increasing cou- 
pling parameter for frequencies up to around 500 cm . 
This is in qualitative agreement with experiment. For 
higher frequencies, the model calculations predict a sharp 
decrease in the coupling parameter. This is not confirmed 
by experiment. 
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The quantitative agreement with experiment is rather 
poor, especially in the low-frequency regime; both mod- 
els provide substantially less activity in this regime than 
observed in experiment. For this reason, Raman spectra 
calculated using either of the two models yield TA peak 
amplitudes far below experimental values. This, in our 
opinion, is their major shortcoming. 

This point is further illustrated in the top graph of 
Fig. 2, where we show the Raman spectrum calculated 
using the MZ model (solid) and an experimental spec- 
trum (dashed) taken from Ref. [21]. The experimental 
spectrum was obtained from ion-implanted o-Si which 
had been annealed at 500°C for two hours. Agreement 
between model and experiment, particularly in the low- 
frequency regime, is poor. 
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coupling parameter [21,22]. The results of this approach 
are illustrated in the bottom graph of Fig. 2. Here, we 
show the Raman spectrum obtained by multiplying a 
computer-generated VDOS with experimental coupling 
parameter data (solid). To calculate the VDOS, we used 
a well-annealed a-Si sample with AO = 10.0°; the experi- 
mental coupling parameter was taken from Ref. [16]. The 
dashed line shows again the experimental Raman spec- 
trum taken from Ref. [21]. Agreement with experiment 
has improved substantially. 

Of the two models considered here, the MZ model pro- 
vides slightly more activity in the low-frequency regime 
than the model of Alben; comparison with experiment 
would therefore favor the MZ model. However, given the 
overall poor performance of both models, we will also 
show in the following sections results obtained using the 
semi-experimental approach. 



B. Raman TO peak width vs. AO 

The intensity of the Raman TO peak decreases 
abruptly on the high frequency side, but not on the low 
frequency side. Beeman therefore defines T as twice the 
half-width at half the maximum height on the high fre- 
quency side of the TO peak, as a meaningful parameter 
to specify the TO peak width [6]. In this paper, we use 
the same definition. 

The solid lines in Fig. 3 show the relation between T 
and A8 for HV polarized light, derived using the model 
of Alben (top), the MZ model (middle) and the semi- 
experimental approach (bottom). Also shown is the re- 
sult obtained by Beeman (dashed). 
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FIG. 2. top: Reduced HV Raman spectrum for a-Si 
calculated using the MZ model (solid), bottom: Reduced 
HV Raman spectrum calculated using the semi-experimental 
approach; the VDOS is obtained by simulation, the cou- 
pling parameter is taken from Ref. [16]. The dashed line in 
both graphs shows an experimental reduced Raman spectrum 
taken from Ref. [21]. Frequency is in cm -1 and all spectra 
are area-normalized. 

Given the overall poor performance of both the model 
of Alben and the MZ model, it may be feasible to fol- 
low a semi-experimental approach in which a computer 
generated a-Si sample is used to calculate the VDOS and 
experimental data is used to describe the coupling pa- 
rameter. This approach is justified because several stud- 
ies have indicated that changes in the Raman spectrum 
are due to changes in the VDOS and not to changes in the 
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FIG. 3. HV Raman TO peak width T/2 as a function of 
AO, calculated using the model of Alben (top), the MZ model 
(middle) and the semi-experimental approach (bottom) . The 
solid lines are linear least-squares fits; the dashed line is the 
result of Beeman: T/2 = 7.5 + 3A0. The units of T/2 and AO 
are cm -1 and degrees, respectively. 

Both the model of Alben and the MZ model produce 
similar results; linear least square fits yield the equations 
T/2 = 3.0A6> - 6.5 and T/2 = 3.7 AO - 13.3, respectively. 
Compared to the result of Beeman, T/2 = 3.0A0 + 7.5, 
we see agreement on the sensitivity (i.e. slope of the lines) 
of T to A9, but not on the overall offset (i.e. intercepts of 
the lines). The same holds for the result obtained in the 
semi-experimental approach; least square fitting yields 
T/2 = 3.3A0 + 9.2 in that case. 



C. Raman TO peak position vs. AO 

As another way to obtain structural information on a- 
Si from its Raman spectrum, we investigate the relation 
between the TO peak frequency (ojto) and AO. Fig. 4 



shows the relation between u>to an d AO for HV polar- 
ized light, derived using the model of Alben (top), the 
MZ model (middle) and the semi-experimental approach 
(bottom). 



505 

500 
495 
490 
485 



480 





HV Alben 



























10.0 11.0 12.0 13.0 14.0 
A6 



505 
500 
495 
g 490 
485 
480 
475 





HVMZ 








































470 



465 



10.0 11.0 12.0 13.0 14.0 
A6 



FIG. 4. Reduced HV Raman TO peak position as a func- 
tion of AO for o-Si, calculated using the model of Alben (top), 
the MZ model (middle) and the semi-experimental approach 
(bottom). The peak position and AO are given in cm -1 and 
degrees, respectively. The solid lines are least-squares fits 
through the datapoints. 

According to Fig. 4, loto shifts to higher frequency 
as AO decreases. However, agreement with experiment 
for both the model of Alben and the MZ model is poor. 
Experimental Raman spectra of well-annealed a-Si sam- 
ples yield ujto m the order of 480 cm -1 [21]. For well- 
annealed configurations, for which AO ranges from 9.9 to 
11.0 degrees, the models exceed the experimental value 
by around 20 cm -1 . 

The semi-experimental approach is in much better 
agreement with experiment; a linear fit yields the equa- 
tion uj to = -2.5A0 + 505.5 which for AO = 10.0° leads 
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to ujto — 480.5 cm -1 . 

The cause of this is that in the model of Alben and the 
MZ model, ojto is determined by the peak in the coupling 
parameter, whereas in the semi-experimental approach, 
it is determined by the VDOS. 



D. Raman TA/TO intensity ratio vs. AO 

Next, we confirm that the TA/TO intensity ratio is 
directly related to AO. Fig. 5 shows the relation between 
reduced Raman TA/TO intensity and AO for HV polar- 
ized light for the model of Alben (top), the MZ model 
(middle) and the semi-experimental approach (bottom) 

From Fig. 5, we see that the TA/TO intensity ratio 
increases with increasing Ad. The increase is approxi- 
mately linear. 
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FIG. 5. HV Raman TA/TO intensity ratio as a function 
of AO for a-Si, calculated using the model of Alben (top), 
the MZ model (middle) and the semi-experimental approach 
(bottom). The solid lines are least-squares fits. AO is in 
degrees. 



HV Raman experiments on well-annealed a-Si samples 
yield a reduced TA/TO intensity ratio around 0.11 [21]. 
For both the model of Alben and the MZ model, the HV 
TA/TO ratio is of order 10 -2 , i.e. one order of magni- 
tude below the experimental value. This is consistent 
with the earlier finding that both models underestimate 
the Raman activity in the low-frequency regime of the 
spectrum, see section III A. 

The semi-experimental approach, on the other hand, 
yields a TA/TO ratio of 0.14 for A9 = 10.0°, in better 
agreement with experiment. 



IV. CONCLUSIONS 

We have generated nine independent sets of 1000-atom 
samples of a-Si that display a variety of short-range or- 
der; the spread in nearest-neighbor bond angles ranges 
from 10 to 14 degrees. For these samples, the HV Raman 
spectra are calculated. To describe the Raman scatter- 
ing process, we have used the earlier bond polarizibility 
model of Alben et al., the more recent model of Mari- 
nov and Zotov as well as experimental data taken from 
Ref. [16]. 

Comparison to experiment shows that both the model 
of Alben and the MZ model greatly underestimate the 
Raman activity in the low-frequency regime of the spec- 
trum. This makes these models less suitable to describe 
low-frequency features of the Raman spectrum, for in- 
stance the TA peak. Of the two models considered here, 
the MZ model is closer to experiment. However, for a 
more accurate calculation of Raman spectra, we propose 
a semi-experimental approach. In this approach, the 
VDOS is obtained in computer simulation and experi- 
mental data is used to describe the coupling parameter. 

As ways to obtain structural information on a-Si from 
its Raman spectrum, we have investigated the relation 
between the TO peak-width T and A0 1 as well as the 
relations between the TO peak frequency and the TA/TO 
intensity ratio as functions of AO. 

According to our results, where we used the semi- 
experimental approach, T and A9 are related by T/2 = 
3.3A0+9.2 for HV polarized light. Here, T is in cm -1 and 
AO in degrees. Comparing this to the result of Beeman 
(r/2 = 3A# + 7.5), we find that our result is similar. 

Our results also show a shift of the Raman TO peak 
frequency (loto) towards higher frequency, as A9 de- 
creases. In the semi-experimental approach, a linear 
least-square fit yields loto — — 2.5A0 + 505.5 for HV 
polarized light. Here, loto is in cm -1 and AO in degrees. 
According to this equation, the shift of ujto is approxi- 
mately 7.5 cm^ 1 , going from unannealed a-Si {AO ~ 13°) 
to annealed a-Si {AO w 10°). This is in quantitative 
agreement with experiment [21]. 

Finally, we have shown that the reduced Raman 
TA/TO intensity ratio (7) is directly related to AO; 
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I decreases linearly with decreasing Ad. Using the 
semi-experimental approach, we obtain the relation / = 
0.0078A6* + 0.0606, where A6 is in degrees. 
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